Long-range correlations and ensemble inequivalence in a generalized ABC model 
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A generalization of the ABC model, a one-dimensional model of a driven system of three particle 
species with local dynamics, is introduced, in which the model evolves under either (i) density- 
conserving or (ii) nonconserving dynamics. For equal average densities of the three species, both 
dynamical models are demonstrated to exhibit detailed balance with respect to a Hamiltonian with 
long-range interactions. The model is found to exhibit two distinct phase diagrams, corresponding 
to the canonical (density-conserving) and grand canonical (density nonconserving) ensembles, as 
expected in long-range interacting systems. The implication of this result to nonequilibrium steady 
states, such as those of the ABC model with unequal average densities, are briefly discussed. 
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In many cases nonequilibrium systems driven by an ex- 
ternal field, such as the temperature gradient or electric 
field, reach a nonequilibrium steadystate which is charac- 
terized by long-range correlations flj . In some cases such 
long-range correlations lead to phenomena like the emer- 
gence of long-range order and spontaneous symmetry 
breaking, even when the dynamics is local and stochas- 
tic 0. This is particularly intriguing in one-dimensional 
models, where such phenomena do not exist under equi- 
librium conditions for systems with short-range interac- 
tions at finite temperatures. 

A model which is often discussed in this context is the 
ABC model This is a one-dimensional three species 

model of a driven system with local dynamics which ex- 
hibits phase separation at nonvanishing densities of the 
three species. It has been shown that at a particular 
point of its phase space, namely, for equal densities of 
the three species, the model obeys detailed balance, and 
the steady-state distribution can be expressed in terms 
of an effective Hamiltonian with long-range interactions. 
Thus, although the dynamics is local and stochastic, the 
long-range nature of the resulting effective interactions 
leads to phase separation and long-range order. 

Equilibrium systems with short-range interactions typ- 
ically evolve into a unique state which is independent of 
their dynamics. Thus, for instance, the Ising models with 
magnetization-conserving (Kawasaki) dynamics and non- 
conserving (Glauber) dynamics result in the same equi- 
librium state. This equivalence of ensembles is broken 
in systems with long-range interactions, where the two- 
body potential decays at large distance with a power 
smaller than the spatial dimension. Such systems are 
nonadditive and their thermodynamic functions are su- 
perextensive with respect to their size. As a result, they 
exhibit unusual phenomena such as ensemble inequiva- 
lence, negative specific heat, and slow relaxation dynam- 
ics (see, for example, The ABC model is expected 
to share some of these properties at the equal densities 
point, where it behaves as an effective long-range inter- 
acting system. A different mechanism for ensemble in- 



equivalence has been discussed within the context of the 
zero-range process, whereby long-range order emerges in 
one dimension due to the noncompactness of the local 
order parameter [9(. 

In this Letter, we generalize the ABC model to include 
vacancies and nonconserving processes. This allows us to 
compare the phase diagram of the conserving model with 
that of the nonconserving model. We show that for equal 
densities both dynamics obey detailed balance with re- 
spect to a Hamiltonian with long-range interactions. Un- 
der conserving dynamics, where the nonconserving pro- 
cesses are excluded, the resulting steady state is that of 
the canonical ensemble. The nonconserving dynamics 
correspond to the grand canonical ensemble. Evaluat- 
ing the phase diagram of the generalized model, we find 
that, as is common in systems with long-range interac- 
tions, the canonical and the grand canonical ensembles 
are inequivalent, yielding different phase diagrams. 

Small deviations from the equal densities condition are 
not expected to significantly alter the steady state. Thus, 
a detailed study of the model with equal densities may 
serve as a guideline for investigation of the general case 
with nonequal densities, where detailed balance is not 
satisfied and an effective energy cannot be defined. This 
may shed some light on the mechanisms leading to long- 
range phenomena in nonequilibrium steady states. 

We begin by outlining the main features of the ABC 
model on a ring [3]. The model consists of three species 
of particles, labeled A, B, and C , which occupy the sites 
of a periodic lattice of length L. The number of particles 
of each type is given by Na, Nb, and Nq with Na + 
Nb + Nq = L. The model evolves by local, random 
sequential dynamics whereby two neighboring particles 
are exchanged clockwise with the following rates: 



AB <=* BA ; BC <=* CB ; CA AC. 
ill 



(1) 



For q = 1 the dynamics is symmetric and the sys- 
tem relaxes to a homogeneous state, in which all par- 
ticles are uniformly distributed, regardless of type. For 
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q 1 the system reaches a nonequilibrium steady state 
in which the particles phase separate into three do- 
mains. The domains are arranged, clockwise, in the order 
AA . . . ABB . . . BCC . . . C for q < 1 and counterclock- 
wise for q > 1. We will assume q < 1 for the rest of 
the Letter. It has been shown that at equal densities, 
Na = Nb = Nc = L/3, the dynamics obeys detailed 
balance with respect to the Hamiltonian 



L-l L- 



i+jJ 



(2) 



= 1 3=1 



Here {JQ} = {^Bt, Q}, 



. _ f 1 if site i is occupied by an A particle , , 
1 y otherwise, 

and similarly for Bi and Cj. The change in H due to the 
exchange of any pair of neighboring particles as described 
in (TTJ) is AH = ±1, where, for example, the exchange of 
AB to BA raises the energy by 1. Detailed balance is 
thus maintained with respect to the steady-state distri- 
bution P({Xi}) = q n ^ x '^ /Z Ll with the partition sum 
Zl = Yl{x-\ q^^ X 4 • While the microscopic dynamics 
((T|) is strictly local, the interactions in T-L are long-ranged 
(mean- field like). 

It has been demonstrated that, by considering an In- 
dependent q which approaches 1 sufficiently fast at large 
L, the q = 1 homogeneous state may be realized [4(. 
For q = exp (— /3/L), the model exhibits a phase tran- 
sition from a homogeneous state for j3 < f3 c to a phase- 
separated state for (3 > f3 c , with /3 C — 2nV3 . The pa- 
rameter /3 serves as an inverse temperature. 

We now generalize the model, allowing us to compare 
conserving and nonconserving dynamics. As a first step 
we introduce vacancies into the model. Thus, each site 
may be occupied either by one of the three species or by 
a vacancy 0. Hence, N = N A + N B + N c < L. The 
dynamics is such that vacancies are neutral, so that a 
particle of any species may hop to the left or right into 
an empty site with equal probability. The following rule 
is added to the exchange rules Eq. (Q}: 



X0 OX , X = A, B, C. 
i 



(4) 



Next, we add nonconserving processes, whereby parti- 
cles are allowed to leave and enter the system in ordered 
groups of three neighboring particles: 



ABC 



000, 



(5) 



where p is a chemical potential, taken to be equal for 
all three species, and p is a parameter whose value does 
not affect the steady state in the case where detailed bal- 
ance is satisfied. This particular choice of the noncon- 
serving process maintains equal densities whenever the 
initial configuration satisfies this condition. 



Focusing on the case of equal densities, we consider 
two alternative dynamics: (i) density-conserving dynam- 
ics where the evolution takes place by the processes ([T]) 
and @ and (ii) nonconserving dynamics where all pro- 
cesses (TTJ, ([J| and §5§ are allowed. Under both dynamics, 
detailed balance is satisfied with respect to the following 
Hamiltonian: 



n ({x 4 }) = u ({Xi}) 



N(N-l) 
6 



fj.NL. (6) 



In the conserving case, the chemical potential term is 
constant and may be omitted. Note that, while the gen- 
eralized ABC dynamics is strictly local, the Hamiltonian 
is long-ranged. The fact that detailed balance is sat- 
isfied with respect to the particle-conserving processes 
(TTJ) and (0J is a result of the neutrality of the vacan- 
cies. The reason why detailed balance is also satisfied 
with respect to the nonconserving processes §5§ has to 
do with the fact that the energy E of a configuration is 
invariant under any translation of ABC triplets, namely, 
E(. . . Y ABC . . .) = E{. . . ABCY . . .), where Y stands 
for either a particle of any species or a vacancy. Thus, 
the rates of depositing or evaporating ABC triplets could 
be taken as independent of the microscopic configuration 
in which these processes take place, as given by the rates 
([S]) . The change in energy corresponding to depositing an 
ABC triplet in a system containing N particles is given 
by Aii = N + 1 - \ (6JV + 6) - 3fiL = -3/xL. Thus, 
detailed balance is satisfied with respect to ([S])- 

The conserving dynamics leads to the canonical steady 
state of the Hamiltonian ([B]), while the nonconserving 
dynamics leads to the grand canonical one. As a re- 
sult of the long-range nature of the interactions, the 
two ensembles need not be equivalent, yielding different 
phase diagrams. In order to carry out this analysis we 
take q = exp(— /3/L) and turn to the continuum limit 
The local concentration of A, B, and C particles 
at the point x = i/L is represented by the density profile 
p n {x) (n = A,B, C) with p(x) — pa{x) + pb{x) + pc (x) ■ 
The steady-state distribution of the density profiles is 
given by P [p n {x)} = expl-LJ 7 [p n (x)}}, where T [p n {x)] 
is the free energy functional. The equilibrium profile can 
thus be found by minimizing the free energy functional 
with respect to p n (x) under the equal densities condition 
/ p n (x)dx = r/3. 

We begin by considering the conserving dynamics (i), 
where the overall density r is fixed, and the free energy 
functional is 



F[pn{x) 



dx[p a(x) hi p A (x) + pb{x) lnp B (x) 



+pc{x)lnp c {x) + (l-p(aO)ln(l-/>(a;))] 
+/3 dx J dz[p A (x)pc(x + z) 

+p B (x)p A (x + z) + p c (x)p B (x + z)] - V) . (7) 
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Here the first term corresponds to the entropy, and the 
second term corresponds to the energy. 

At high temperatures T = 1//3, the free energy is min- 
imized by the homogeneous density profiles p n (x) = r/3, 
corresponding to the disordered phase. In order to find 
the transition, we note that in the ordered phase the den- 
sity profiles are expected to satisfy pb{x) = pa(x — 1/3) 
and pc = Pa(x + 1/3). Assuming a smooth transition to 
the ordered phase we expand pa(x) close to the homoge- 
neous solution: 



Pa(x) 



a m cos(27rma;) 



(8) 



This profile is arbitrarily chosen among all its transla- 
tionally related degenerate profiles. Substituting ([5]) in 
Eq. ([7]), together with the symmetry conditions on Pb(x) 
and pc(x), results in a series expansion for the free en- 
ergy J- in powers of the amplitudes a m . Instabilities of 
the homogeneous state are dominated by the first mode 
Oi, while higher modes driven by it. Applying the 

equilibrium conditions dJ 7 /da rn — 0, one can express the 
coefficients a m for m > 1 in powers of a%. This leads to a 
Landau expansion for the free energy J 7 , with a\ serving 
as an order parameter: 



F[ PA {x)\ = T (~) + ha\ + U4 + 



where 



/a 
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4r 



3^/3 
8tt 



81 
32^ 



6ir 



12tt 



(9) 



(10) 



Since f± > for any r and /3, there is a second-order 
phase transition at fa = 0, where j3 c — 2i:^/?>/r. 

By evaluating the chemical potential using p = 
jj Jj^J 7 (r / 3) , the critical line is obtained: 



1 



In 1 



(11) 



The (l//3,/x) phase diagram is shown in Fig. [TJ 

In order to evaluate the phase diagram correspond- 
ing to the nonconserving dynamics (ii), we consider the 
Gibbs free energy Q[p n (x)] — T \p n {x)\ — j3pr and study 
the stability of the homogeneous phase at a given chemi- 
cal potential p. To this end we expand Q in small devia- 
tions from the homogeneous profile. Here, however, one 
should also allow for variations of the overall density Sr. 
Thus, the A-particle density profile close to the transition 
can be written as 

S x 

Pa{x) = g + y + ^2 a m cos(2irmx). (12) 

m— 1 

As before, the B and C density profiles are obtained 
by translation operations. The equilibrium conditions 
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Figure 1: The (1//3, p) phase diagram of the ABC model. For 
conserving dynamics the model exhibits a line of continuous 
transition [thin solid line, Eq. {TTJ]. For nonconserving dy- 
namics the model exhibits a multicritical (fourth-order) point 
(star), beyond which the transition becomes first-order (thick 
solid line) determined numerically (circles). 



dQ jd (5r) — and dQ /da m = 0, for m > 1 result in the 
Landau expansion of Q in terms of oi, 

g [p A {x)\ = g (0 + 9 2a\ + gi a\ + 9& a\ + ... , (13) 



with 



.92 = Ii 



3V3/3 



4r 



(14) 



This coefficient vanishes at /3 C = 27rv / 3/?~, yielding the 
critical line, as long as 54 > 0. On the critical line the 
fourth-order coefficient is given by 



94 {Pc) 



27 
32^3 



(3r-l) 



(15) 



The transition is, thus, continuous for r > 1/3, becoming 
first-order at a multicritical point (MCP) r = 1/3, where 
.92 = .94 = 0. In the (l/f3,p) plane, the MCP is given by 



/3mcp = GttVs ; pmcp = - 



In 6 
6ttV^ 



-0.0549. (16) 



Calculating higher-order terms in the expansion f|13f) , we 
find that at the MCP g§ = while g$ > 0, which implies 
that this is a fourth-order critical point. 

To complete the phase diagram we note that at T = 
the system is fully phase-separated, with randomly dis- 
tributed vacancies. For such states one has Gt=o — 
—r 2 /18 — pr. Minimizing Q with respect to r, one finds a 
first-order transition at p = —1/18 between an empty 
state (r = 0) for p < —1/18 and a phase-separated 
state with r = 1 for p > —1/18. The first-order line 
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Figure 2: The r(fj,) curve obtained from Monte Carlo simu- 
lations with L = 1800 and temperatures (a) above and (b) 
below the multicritical point. In (a) the two types of dynam- 
ics result in the same curve with a continuous transition. The 
expected critical point is marked by x; deviations are due to 
the finite size of the system. In (b) the conserving dynam- 
ics results in a continuous transition (x), while the noncon- 
serving dynamics displays a discontinuity in r at a first-order 
transition. Hysteretic behavior is observed. The expected 
transition point obtained by minimization of Q is indicated 
by the vertical line. Dotted lines correspond to homogeneous 
profiles. 



connecting the MCP and the T — transition point is 
found numerically by integrating the dynamical equa- 
tions dp n (x,t)/dt = —5Q/Sp n (x,t). While these equa- 
tions do not represent the actual dynamics of the system, 
their steady-state solution reproduces the steady-state 
profile corresponding to the minimum of Q. 

The phase diagram of the nonconserving ABC model 
is given in Fig. [TJ where it is compared with that 
of the conserving model. While the transition to the 
phase-separated state in the conserving model is con- 
tinuous throughout the (1//3, /i) plane, the transition in 
the nonconserving model changes character at a fourth- 
order critical point. The transition lines in the two mod- 
els coincide when both are second-order. However, the 
nonconserving model tends to enhance phase separation, 
where it is found to exist in domains in the fi) plane 
in which the conserving model displays a homogeneous 
phase. This is a characteristic behavior of systems with 
long-range interactions, resulting from the nonadditivity 



of these systems @, H| ■ 

In order to verify the picture emerging from the contin- 
uum approximation, Monte Carlo simulations were per- 
formed for (i) conserving dynamics [rules ([lj and (QJ] and 
(ii) nonconserving dynamics [rules ([IJ, (0| and fl5J], at 
temperatures above and below the MCP (Tmcp — 0.03). 
In the density-conserving case, the chemical potential 
fi is calculated by using the Creutz algorithm In 
this method, particles may be exchanged between the 
system and an external single degree of freedom, a "de- 
mon," with Nsystem + N demon = N. The chemical po- 
tential /i is determined by calculating the distribution: 
P(N demon ) ~ exp (-/3^N demon ). The results for the 
density curves r(/t) are shown in Fig. [5] Above the 
MCP, both the conserving and nonconserving dynamics 
yield the same curve within the numerical accuracy. On 
the other hand, below Tmcp the nonconserving dynamics 
leads to a discontinuous change in the total density, with 
hysteretic behavior, consistent with a first order phase 
transition. The conserving dynamics result in a contin- 
uous transition at the expected critical point, realizing 
density values which cannot be reached under noncon- 
serving dynamics. 

The canonical and grand canonical phase diagrams of 
the generalized ABC model at equal densities may serve 
as a very useful starting point for analyzing the phase 
diagram at nonequal densities, where detailed balance is 
not satisfied and where a free energy cannot be defined. 
Indeed, analysis of the dynamical equations in the contin- 
uum limit and numerical simulations carried out with the 
density of one species different from the other two yield 
a qualitatively similar phase diagram to that of Fig. [1] 
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